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Abstract 

This  paper  shows  how  /i-node,  e-edge  graphs  can  be  contracted  in  a  manner  similar  to  the 
parallel  tree  contraction  algorithm  due  to  Miller  and  Reif.  We  give  an  0((n  +  e)/lg  n)- 
processor  deterministic  algorithm  that  contracts  a  graph  in  0(lg^)  time  in  the  EREW 
PRAM  model.  We  also  give  an  0(/i/lg  n)-processor  randomized  ^gorithm  that  with  high 
probability  can  contract  a  bounded-degree  graph  in  0(lg  n  +  Ig^  7)  time,  where  7  is  the 
maximum  genus  of  any  connected  component  of  the  graph.  (The  algorithm  can  be  made  to 
run  in  deterministic  0(lg  n  \gn  +  lg27)  time  using  known  techniques.)  This  algorithm  does 
not  require  a  priori  knowledge  of  the  genus  of  the  graph  to  be  contracted.  The  contraction 
algorithm  for  bounded-degree  graphs  can  be  used  directly  to  solve  the  problem  of  region 
labeling  in  vision  systems,  i.e.,  determining  the  connected  components  of  bounded-degree 
planar  graphs  in  0(lg  n)  time,  thus  improving  the  best  previous  bound  of  0(lg2  n). 
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Abstract 

This  paper  shows  how  n-aode,  e-edge  graphs  can  be  eon- 
traeted  in  a  manner  similar  to  the  parallel  tree  contrac¬ 
tion  algorithm  due  to  Miller  and  Reif.  We  give  an  0((n  -f 
e)/lgn)-processor  deterministic  algorithm  that  contracts 
a  graph  in  0(lg^  n)  time  in  the  EREW  PRAM  model. 
We  also  give  an  0(n/  Ig  n)-proce88or  randomized  algorithm 
that  with  high  probability  can  contract  a  bonnded-degree 
graph  in  0(lg  n+lg^  t)  time,  where  i  is  the  maximum  genus 
of  any  connected  component  of  the  graph.  (The  algorithm 
can  be  made  to  run  in  deterministic  0(lg  n  Ig*  n  +  Ig^  7) 
time  using  known  techniques.)  This  algorithm  does  cot 
require  a  priori  knowledge  of  the  genus  of  the  graph  to 
be  contracted.  The  contraction  algorithm  for  bonnded- 
degree  graphs  can  be  used  directly  to  solve  the  problem 
of  region  labeling  in  vision  systems,  i.e.,  determining  the 
connected  components  of  bounded-degree  planar  graphs  in 
O(lgn)  time,  thus  improving  the  best  previous  bound  of 
0(lg»n). 


1  Introduction 

The  parallel  tree  contraction  technique  of  Miller  and  Reif 
[2^  has  proved  to  be  a  valuable  tool  in  many  parallel  graph 
algorithms.  This  paper  provides  a  similar  contraction  tech¬ 
nique  that  applies  not  only  to  trees,  but  also  to  general 
graphs.  It  also  provides  a  second,  potentially  faster  con¬ 
traction  technique  for  bounded-degree  graphs  (and  general 
graphs  when  an  embedding  is  known). - - - — 

A  contraction  $tep  of  an  undirected  graph  G  s  {V,E) 
with  respect  to  an  edge  (v.v)  £  E  u  the  operation  that 
replaces  «  and  o  by  a  new  vertex  w  which  is  adjacent  to  all 
those  vertices  to  which  •  and  v  were  adjacent.  If  vertices  m 
and  V  have  an  adjacent  vertex  x  in  common,  the  contrac¬ 
tion  step  produces  only  one  edge  between  w  and  x  rather 
than  two.  Thus,  the  graph  that  results  from  a  contraction 
step  is  not  a  multigraph.  A  parallel  paired  contraction  etep 

^'This  research  is  supported  in  part  by  the  Delame  Advanced 
Reaeardb  Projects  Agency  under  Contract  N00014-87-K-083S 
and  in  part  by  the  Office  of  Naval  Rwsesrcb  under  Contract 
N00014-86-K-0593  and  in  part  by  an  IBM  graduate  fellowehip. 


is  a  vertex-independent  sequence  of  contraction  steps.  A 
parallel  malti-eontraetion  itep  is  a  sequence  of  contraction 
steps  with  respect  to  an  acyclic  set  of  edges.  We  refer  to  the 
sequence  as  simply  a  parallel  contraction  step  whenever  the 
type  of  contraction  is  clear  bom  context.  A  parallel  paired 
contraction  step  of  any  bonnded-degree  graph  can  be  imple¬ 
mented  in  constant  parallel  time.  A  (par^lel)  contraction 
of  a  graph  is  the  process  by  which  a  connected  graph  is 
reduced  to  a  single  node  by  iterated  (parallel)  contraction 
steps. 

Miller  and  Reif  [23]  show  how  any  tree  can  be  contracted 
in  0(lg  n)  time  nsing  0(n)  processors  in  the  CRCW  PRAM 
model  using  randomization.  (They  also  show  how  to  re¬ 
duce  the  number  of  processors  to  0(n/lgn))  They  give  an 
0(lgn)-time  algorithm  for  tree  contraction  that  is  deter¬ 
ministic,  but  the  deterministic  algorithm  does  not  perform 
contraction  in  the  strict  sense  described  above.  Tree  con¬ 
traction  can  be  made  to  run  in  randomized  O(lgtt)  time 
and  in  deterministic  0(lg  n  Ig*  n)  time  using  0(n)  proces¬ 
sors  in  the  EREW  and  DRAM  models,  as  was  shown  in 
(20). 

In  this  paper/  we  use  the  technique  of  parallel  contrac¬ 
tion  in  a  more  general  setting.  We  show  that  a  remarkably 
ample  contra^on  algorithm  that  uses  (n  -f  e)/  Ig  n  proces¬ 
sors  rednces  k  connected  n-node  e-edge  graph  to  a  single 
node  in  0(lg*  w)  steps  and  a  second  sim^e  algorithm  which 
naes  the  first  as  a  subrouting  red  aces  a  connected  bound ed- 
degree  graph  to  a  single  node  using  w/  Ig  w  pTSceisors  in 
0(lgn-f  Ig’  7)  steps,  where  7  is  the  maximum  genus  of  any 
connected  component  of  graph  G.  The  second  algorithm 
immediately  yidds  an  asymptotically  efficient  solution  to 
the  problem  of  region  labeling  in  vision  systems,  as  well  as 
to  solutions  of  other  planar  graph  problems. 

The  genus  of  a  graph  is  the  sum  of  the  geni  of  the  con¬ 
nected  components.  That  is,  if  the  tth  connected  com- 
TOxent  has  genus  7.,  then  the  genus  of  the  graph  is  70  = 
V .  7i.  The  running  time  of  our  second  algorithm,  however, 
Spends  only  upon  the  maximum  genus  of  any  connected 
component,  (i.  e.  7  s  maxi7i).  Henceforth,  when  we  refer 
to  the  genus  of  a  graph  G,  we  will  use  the  notation  7a  when 
we  mean  the  true  genus  and  we  will  use  the  notation  7  when 
we  mean  the  maximum  genus  of  any  connected  component 
of  graph  G. 

Figure  1  illustrates  the  data  structure  we  use  in  our  algo¬ 
rithms.  We  represent  each  vertex  of  degree  d  by  a  doubly- 
linked  ring  of  d  real  processors  alternating  with  d  dummy 
processors.  The  dummy  processors  are  needed  for  efficient 
imidementation  of  a  parallel  multi-contraction  step.  The 
edges  between  real  processors  and  dummy  processors,  rep-  i 
resented  by  heavy  lines  in  figure  1 ,  are  called  vertex  edges 
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Figure  1:  We  represent  graphs  using  a  ring  of  proces¬ 
sors  for  each  vertex.  Two  (round)  real  processors  are 
separated  by  a  (square)  dummy  processor  in  the  list. 
Vertex  edges,  drawn  as  heavy  lines,  go  between  proces¬ 
sors  in  the  same  vertex  ring,  while  graph  edges  (thin 
lines)  go  between  processors  in  different  vertex  rings. 

and  edges  between  vertex  rings  are  called  graph  edgei.  Each 
processor  in  a  vertex  ring  knows  the  ID  of  the  vertex,  de¬ 
fined  to  be  the  maxiniain  identifier  of  any  processor  in  the 
vertex  ring. 

The  remainder  of  the  paper  is  organized  as  follows. 
Section  2  presents  the  contraction  algorithm  for  general 
graphs,  argues  its  correctness,  and  analyses  its  running 
time.  Section  3  presents  the  contraction  algorithm  for 
bounded-degree  graphs  and  Section  4  analyzes  its  running 
time  using  a  'hnissing  edge”  lemma  which  is  proved  in  sec¬ 
tion  5.  Section  6  shows  how  graph  contraction  can  be 
applied  to  various  graph  problems,  including  the  region¬ 
labeling  problem.  Section  7  offers  some  concluding  re¬ 
marks. 

2  General  Algorithm 

This  section  presents  an  0(n  -f  e)-processor  contraction  al¬ 
gorithm  for  general  graphs,  argues  its  correctness  and  an¬ 
alyzes  its  running  time.  We  describe  how  to  contract  an 
edge  and  we  give  a  strategy  for  determining  a  set  of  acyclic 
edges  to  contract  in  a  parallel  multi-contraction  step.  We 
then  discuss  implementation  issues  such  as  avoiding  sort¬ 
ing,  maintaining  linear  space,  and  reducing  the  processor 
count  to  0((n  +  e)/  Ig  n). 

Figure  2  illustrates  how  to  contract  an  edge  («,  v)  by 
using  the  dummy  processors  to  merge  the  two  vertex  rings 
and  splice  out  the  real  processors  corresponding  to  the  con¬ 
tracting  edge.  The  dummy  edges  allow  us  to  simultane- 


Figure  2;  To  contract  a  graph  edge  (u,  v),  we  use  the 
neighboring  dummy  processors  on  each  side  to  merge 
the  vertex  rings  and  splice  out  the  real  processors.  We 
later  clean  up  the  extra  dummy  processors.  Dummy 
processors  allow  contraction  of  trees  of  arbitrary  depth. 

ously  contract  any  acyclic  set  of  edges,  including  trees  of 
arbitrary  depth,  in  constant  time.  We  then  pay  an  extra 
O(lgn)  time  to  remove  the  extra  dummy  vertices. 

The  algorithm  for  contracting  graphs  is  simple.  It  uses  a 
subroutine  MERaE(n,o)  which  performs  the  ^ge  contrac¬ 
tion  described  above  for  edge  (n, «). 

CONTJUCT 

1  While  3  a  vertex  with  deg  ee  >  0  do  In  parallel 

2  Each  vertex  «  with  degree  >  0,  do  in  parallel 

3  Let  V  be  the  neighbor  of  highest  ID 

4  Meroe(u,  v) 

5  Remove  extra  dummy  processors 

6  Remove  multiple  edges  to  adjacent  vertices 

7  Propagate  new  ID. 


We  now  argue  the  correctness  of  Algorithm  Contract. 
The  only  subtle  point  is  that  simultaneous  contraction  of  a 
set  of  edges  that  form  a  tree  results  in  a  single  vertex  ring. 
If  we  were  to  contract  edges  that  form  a  cycle,  then  we 
may  end  up  with  multiple  rings,  thus  severing  a  connected 
component.  Because  each  vertex  chooses  to  merge  with  its 
nei^bor  of  maximum  ID,  the  edges  chosen  to  contract  in 
an  iteration  cannot  form  a  cycle. 

We  now  analyse  the  running  time  of  algorithm 
Contract.  Each  iteration  of  the  loop  requires  O(lgn) 
time.  The  loop  termination  check  in  line  1  can  be  done 
in  0(lg  n)  time.  Each  processor  can  determine  in  constant 
time  whether  the  vertex  ring  to  which  it  belongs  has  de¬ 
gree  >  0.  We  can  remove  multiple  edges  in  O(lgn)  time 
by  first  sorting  the  processors  in  a  vertex  ring  with  respect 
to  the  ID  of  the  neighbor  across  the  graph  ^ge  and  then 
pointer  jumping.  Similarly,  propagating  the  highest  ID  of 
any  neighbor  to  all  processors  in  a  ring,  removing  extra 
dummy  processors,  and  propagating  the  new  vertex  ID  can 
all  be  done  in  O(lgn)  time  using  pointer  jumping.  If  we 
are  careful,  all  pointer  jumping  can  be  done  without  using 


concurrent  ren^.  The  Merqe  operntios  require*  only  con- 
stnnt  time.  Thu*  each  iteration  require*  0(tg.n)  time.  The 
algorithm  terminates  in  0(lg  n)  iteration*  becanae'each  ver¬ 
tex  merges  on  every  iteration,  at  least  halving  the  nnmber 
of  vertices. 

We  can  avoid  sorting  if  we  allow  multiple  edge*  to  re¬ 
main  after  contraction  and  instead  use  pointer  jumping  to 
remove  the  self  edges  that  appear  in  later  iterations.  If  we 
choose  this  option,  then  in  line  3  we  must  also  break  lies 
among  all  edges  to  vertex  v.  Using  the  idea  of  ‘Spares” 
from  [20],  the  algorithm  can  be  made  to  use  only  linear 
space. 

To  reduce  the  processor  count  to  (n  -f-  e)/lgn,  each 
processor  simulates  k  processors  per  iteration  (initially 
it  =  Ig  n).  After  each  iteration  of  the  contraction  algorithm, 
we  balance  the  work  among  the  processors  by  enumerating 
the  remaining  nodes  in  0(lg  n)  time  via  parcel  prefix  and 
compacting  memory  in  0{k)  time.  Thus  each  iteratioiu 
requires  0(k  Ig  n)  time.  The  parameter  k  is  initially  Ig  n 
and  decreases  exponentially  with  each  iteration,  so  over¬ 
all  the  contraction  algorithm  runs  in  0(lg^  n)  time  using 
(n  -f  e)/lgn  processors.  The  balancing  steps  can  be  done 
without  concurrent  memory  access,  so  these  bounds  hold 
for  the  EREW  PRAM. 


3  Bounded-Degree  Algorithm 

In  this  section  we  present  a  linear-processor  contraction  al¬ 
gorithm  for  bounded-degree  graphs.  We  give  a  randomised 
strategy  for  determining  a  set  of  vertex-disjoint  edges  to 
contract  on  a  given  parallel  contraction  step.  We  discuss 
simplifications  for  the  case  where  the  graph  is  known  to  be 
planar  or  of  low  genus.  We  discuss  implementation  issues 
including  maintaining  linear  space,  using  this  algorithm  for 
arbitrary  graphs  provided  an  embedding  is  known,  mak¬ 
ing  the  algorithm  deterministic,  and  reducing  the  processor 
count.  Finally  we  mention  some  of  the  parallel  models  tor 
which  the  analysis  of  section  4  holds. 

The  contraction  algorithm  for  bounded-degree  graphs  is 
as  follows.  We  assume  that  no  vertex  has  degree  greater 
than  40.  The  number  j  in  line  1  is  a  constant  chosen 
such  that  algorithm  BoundezvContract  terminates  in 
0(lg  n  +  Ig*  7)  time  with  high  probability.  The  constant 
will  be  discussed  during  the  analysis  of  the  algorithm  in 
section  4. 

Algorithm  BoUNDED-CONTRACT 

1  Repeat  j  Ig  n  times 

2  for  each  vertex  u  do  in  parallel 

3  Randomly  choose  an  adjacent  vertex  e  such 

that  MERaE(«, «)  would  produce  a  vertex 
of  degree  at  most  40 

4  for  all  vertices  «  and  v  that  choose  each  other, 

do  Merge(u,v) 

5  for  each  vertex  u  do 

remove  multiple  edges  to  adjacent  vertices 

6  ContractO  ;go  to  general  contraction  algorithm 


The  algorithm  can  be  broken  down  into  two  phases: 
phase  1  in  lines  1-5  in  which  vertices  resulting  from  con¬ 
traction  cannot  have  degree  greater  than  40  and  phase  2  in 
line  6  in  which  vertices  resulting  from  contraction  can  have 
arbitrary  degree.  In  section  4  we  argue  that  phase  1  runs  in 
O(lgn)  time  and  with  high  probability  reduces  an  n- vertex 


genus-7  graph  to  a  graph  whose  largest  connected  com¬ 
ponent  has  0(7)  vertices.  By  the  argument  in  section  2, 
phase  2  reduces  each  0(7)-sued  connected  component  of  a 
graph  to  a  single  vertex  in  0(lg*  7)  time.  Planar  graphs  and 
graph*  of  constant  genus  will  be  reduced  to  constant-sixed 
graphs  during  phase  1  of  algorithm  Bounded-Contract 
with  high  probability.  If  we  know  that  a  graph  is  of  con¬ 
stant  genns,  we  can  simplify  the  algorithm  to  repeated  ap¬ 
plication  of  lines  2-5  only,  nntil  the  graph  is  contracted. 
The  parameter  40  in  line  3  must  be  replaced  by  c(y),  a  con¬ 
stant  depending  upon  the  genus  7.  In  particular,  ^7)  =  40 
suffices  for  planar  graphs.  We  then  check  for  termina¬ 
tion  (all  vertices  degree  0)  every  Ig  n  parallel  contraction 
steps.  This  simplifi^  algorithm  will  contract  an  n-vertex 
bonnded-genus  graph  to  a  single  vertex  in  0(lg  n)  time  with 
high  probabilty.  Fnrthermore,  because  the  degree  of  each 
vertex  remain*  bounded  thronghont  the  algorithm  in  this 
case,  we  can  simply  assign  one  processor  to  each  vertex. 

The  algorithm  can  be  made  to  work  in  0(lgn  +  Ig*  7) 
time  for  any  graph  whose  maximum  degree  is  a  constant 
greater  than  40  by  replacing  the  number  40  in  line  3  with 
the  maximum  degree.  Moreover,  the  algorithm  can  be 
made  to  work  in  0(lg  n  -f  Ig*  7)  time  for  any  graph,  pro¬ 
vided  an  embedding  for  the  graph  is  known,  since  we  can 
transform  any  general  graph  of  genns  7  into  a  degree-3 
genn*-7  graph  by  replacing  each  vertex  by  a  ring  of  degree- 
3  vertices.  Again,  using  the  idea  of  “spares”  from  [20],  the 
algorithm  can  be  made  to  use  only  linear  space. 

The  contraction  algorithm  can  be  modified  to  run  in  de¬ 
terministic  0(lg  n  Ig*  n  -f  Ig*  7)  time  by  nsing  the  0(lg*  n)- 
time  algorithm  of  Goldberg,  Plotkin,  and  Shannon  for  3- 
coloring  rooted  trees  [11]  to  guarentee  that  a  constant 
fraction  of  vertices  merge  on  each  iteration.  Randomness 
is  nsed  only  in  the  choice  of  vertex  pairings  in  phase  1. 
In  the  deterministic  version  of  the  algorithm,  each  vertex 
«  chooses  the  vertex  v  with  smallui  identifier  such  that 
MERaE(u,  u)  would  produce  a  vertex  of  degree  at  most 
dmu  •  the  maximum  degree  of  any  vertex.  As  before,  edges 
chosen  by  both  vertices  adjacent  to  it  are  automatically 
selected  to  contract.  Consider  the  graph  induced  by  the 
edges  chosen  by  exactly  one  adjacent  vertex.  If  the  edges 
are  directed  such  that  the  vertex  choosing  the  edge  is  at 
the  tail,  then  this  graph  is  a  bounded-degree  directed  for¬ 
est  with  edges  directed  from  child  to  parent.  We  can  then 
use  the  Goldberg-Plotldn-Shannon  0(ig*  n)-time  algorithm 
for  coloring  bounded-degree  trees  to  cdor  the  forest  nsing  a 
constant  number  of  color*  [11].  Then,  we  sequence  through 
the  color*  allowing  edges  chosen  by  vertices  of  the  current 
color  to  contract,  provided  the  vertex  on  the  other  end  of 
the  edge  has  not  already  participated  in  a  contraction  dur¬ 
ing  this  parallel  contraction  step.  In  this  scheme,  we  always 
have  a  constant  fraction  of  the  vertices  merging. 

We  can  reduce  the  processor  requirement  of  phase  1  to 
n/  Ig  n  by  nsing  the  techniques  of  Gazit  and  Reif  [10]  who  in 
tarn  use  the  load  balancing  techniques  of  Cole  and  Vishkin 
[5].  Shannon  [26]  gives  a  scheduling  algorithm  that  converts 
the  deterministic  0(lgnlg*  n)-time  n-processor  algorithm 
for  graphs  of  bounded  genns  to  an  optimal  0(lgnlg*n)- 
time  n/(lg  n  Ig*  n)-processor  algorithm. 

.^orithm  Bounded-Contract  is  robust  in  that  it  can 
be  implemented  in  several  of  the  most  restrictive  parallel 
models.  By  careful  attention  to  the  data  structures  used  to 
implement  adjacency  lists,  it  is  possible  to  guarantee  that 
no  concurrent  reading  or  writing  occurs,  and  thus,  the  per¬ 
formance  bounds  apply  in  the  EREW  PRAM  model.  Since 
each  processor  is  responsible  for  a  single  edge  (or  vertex) 


of  the  graph,  U  is  naturally  a  “data-paialiel’’  algorithm  in 
the  sense  of  [16].  Finally,  the  simplified  venion  of  the  al¬ 
gorithm  which  uses  only  phase  1  is  “conserrative”  in  the 
sense  of  [20],  and  thus  runs  in  0(lg  n  +  Ig^  y)  steps  in  the 
DRAM  model. 


traction  step  in  phase  1  of  algorithm  Bounoeo-Contract, 
a  constant  fraction  of  the  vertices  are  eligible  to  contract. 

Lemma  2  Any  gentu-fa  graph  G  =  {V,  E)  with  degree  at 
most  dnu  >  40  and  |V|  >  40^0  has  at  least  ^  |V|  good 
vertices. 


4  Analysis  of  the  Contraction 
Algorithm 

In  this  section  we  analyse  the  contraction  algorithm  of  sec¬ 
tion  3.  We  require  two  lemmas  concerning  bonaded-degree 
graphs.  Lemma  1  (the  Missing- Edge  Lemma)  provides 
an  upper  bound  on  the  number  of  edges  in  a  bounded- 
degree  graph  based  on  Euler’s  formula  and  the  number 
of  degree-three  vertices  that  are  ineligible  for  contraction. 
The  missing-edge  lemma  is  actually  proved  in  section  5. 
We  use  a  pigeonholing  argument  to  show  in  Lemma  2  that 
at  each  parallel  contraction  step  of  phase  1,  a  constant  frac¬ 
tion  of  the  vertices  are  adjacent  to  at  least  one  edge  that 
can  be  contracted.  — This  lemma  is  essentially  the  same  as 
one  proved  independently  by  Boyar  and  Karloff  [3]  in  the 
context  of  coloring  planar  graphs,  but  our  lemma  is  more 
general  and  our  proof  differs  somewhat.  The  final  result  of 
the  section  proves  that  with  high  probability,  the  running 
time  of  Algorithm  Bounded-Contract  is  0(lg  n  -t-  Ig^  t). 

We  analyze  Algorithm  Boundcd-Contract  using  a 
constant  dnu  >  40  in  place  of  40  in  line  3  of  the  con¬ 
traction  algorithm.  Using  a  symbolic  value  allows  ns  to  see 
in  the  analysis  why  we  choose  degree  40  as  the  maximum 
degree  in  the  algorithm  and  to  see  how  the  contraction  algo¬ 
rithm  generalizes  lo  bounded-degree  graphs  in  general.  In 
fact,  a  choice  of  dau  =  3  suffices  for  Algorithm  Bounded- 
Contract  to  contract  binary  trees  in  randomized  O(lgn) 
time. 

We  first  present  some  definitions.  Let  G  =  (V,£)  be 
a  graph  with  degree  at  most  dmu  >  40.  We  call  an 
edge  («,«)  €  E  eligible  if  MERaE(a,t>)  would  produce  a 
vertex  w  of  degree  at  most  daui  as  in  line  3  of  Algo¬ 
rithm  Bounded-Contract.  Typically,  the  degree  of  w 
is  deg(u>)  =  deg(«)  +  deg(o)  —  2,  but  it  is  less  whenever 
the  adjacency  lists  of  s  and  v  have  vertices  in  common  be¬ 
cause  a  contraction  step  removes  multiple  edges.  (In  fact, 
after  the  parallel  contraction  step  implemented  by  Algo¬ 
rithm  Bounded-Contract,  the  degree  of  vertex  tv  may 
be  even  smaller,  since  lines  5  and  6  remove  additional  mul¬ 
tiple  edges  caused  by  other  simultaneous  contraction  steps.) 
We  define  a  vertex  s  €  V  to  be  good  if  it  is  incident  on  at 
least  one  eligible  edge,  and  bad  otherwise.  A  vertex  of  de¬ 
gree  1  or  2  is  automatically  good.  A  vertex  v  of  degree  3 
is  good  unless  it  is  incident  on  three  degree-dnu  vertices 
P],  P2,  and  Vs  that  are  independent  (no  edge  between  any 
pair).  Consequently,  the  bad  degree-3  vertex  v  causes  the 
edges  (oi,vz),  (03,03),  and  (01,03)  to  be  niissinpfrom  the 
graph,  even  though  their  inclusion  would  not  increase  the 
gentu  of  the  graph. 

The  next  lemma  uses  the  notion  of  missing  edges  to  show 
that  a  graph  with  bad  degree-3  vertices  is  sparser  than 
required  by  Euler’s  formula  alone.  It  is  proved  in  section  5. 

Lemma  1  (Missing-Edge  Lemma)  A  graph  G  s  (K,  E)  of 
gentu  fa  with  63  >  0  bad  degree-i  vertices  has  at  most 
SlVj  —  63  IO70  edges. 


Proof.  We  begin  the  proof  by  defining  the  following  sets: 

Vgood  is  the  set  of  good  vertices, 

Vbsds  >•  tbe  set  of  bad  degree-3  vertices, 

Vb^4M  is  the  set  of  bad  vertices  of  degrees  4,  5,  or  6, 
Vtuzk  is  the  set  of  vertices  of  degrees  dm  —  *, 
for  i  s  0, 1,2  or  3, 

=  V  —  (I^od  c  tJ  V^«d456  tJ  I^igh  ). 

having  cardinalities  g,  63,  6434 ,  h,  and  6r«t  respectively. 

We  determine  a  lower  bound  on  the  number  g  of  good 
vertices  using  three  constraints.  The  first  constraint  is  a 
lower  bound  on  the  number  e  of  edges  in  £,  which  is  also 
half  the  sum  of  the  degrees  of  all  vertices  in  V.  For  each 
of  the  sets  defined  above,  we  underestimate  the  sum  of  the 
degrees  of  the  vertices  in  the  set  to  yield 

e  >  —  (g  +  Sis  +  ibiie  (d  max  —  3)  A  + 

=  ^(p  +  36s  +  46«3«  +  (dgux  —  3)6 

7(n  -  p  -  61  -  644#  -  6)).  (1) 

We  can  use  a  similar  technique  to  determine  a  lower  bound 
on  the  number  of  high-degree  vertices.  The  number  of  edges 
leaving  the  set  High  is  at  most  dmu6,  which  must  be  at 
least  the  number  leaving  the  sets  and  Vb»d4ss  since 
each  of  the  vertices  in  these  sets  is  adjacent  only  to  vertices 
in  V^gk.  Hence,  we  obtain 

h  >  +  (2) 

"max 

Finally,  from  the  Missing-Edge  Lemma  we  have 


e  <  3n  —  63  -1-  IOto-  (3) 

We  use  these  three  constraints  to  obtain  the  desired  lower 
bound  on  the  number  g  of  good  vertices.  Combining  In¬ 
equalities  (1)  and  (3)  and  solving  lor  g,  then  substituting 
for  h  using  Inequality  (2)  yields 


>  j(n  -  207a), 


if  dous  >  40  as  assumed.  Since  we  also  assume  that  n  > 
407a,  we  have  that  g  >  n/12  which  completes  the  proof. 


We  now  analyze  the  behavior  of  phase  1.  Phase  1  runs 
in  O(Ign)  time.  Given  a  graph  G  with  maximum  ver¬ 
tex  degree  at  most  40,  Algorithm  Bounded-Contract 
does  not  allow  the  degree  of  any  vertex  to  exceed  40  dur¬ 
ing  phase  1  (the  parallel  contraction  steps  in  lines  2-5  of 
Algorithm  Bounded-Contract).  Therefore,  each  such 
step  can  be  performed  in  constant  time,  since  contraction 
and  removal  of  multiple  edges  involves  only  communicat¬ 
ing  around  constant-length  vertex  rings.  Since  we  execute 
j>  Ign  parallel  contraction  steps  for  a  constant  j,  phase  1  is 
over  in  0(lg  n)  time. 


The  next  lemma  uses  the  Missing- Edge  Lemma  and  pi¬ 
geonholing  to  show  that  during  each  parallel  paired  con¬ 


We  now  use  Lemma  2  to  show  thai  for  euitable  choice 
of  constant  >,  with  high  probability  phaa*  l^ednces  an  n- 
vertex  genas-y  graph  to  a  graph  with  0(7)  vertibea  in  its 
largest  connected  component.  We  need  only  show  that  with 
high  probability,  0(lg  n)  parallel  contraction  steps  suffice 
to  contract  each  connected  component  to  size  0(7).  From 
lemma  2,  we  have  that  at  least  1/12  of  the  vertices  are  good 
provided  that  there  are  at  least  407a  vertices.  The  lemma 
applies  independently  to  each  connected  component.  For 
component  t,  contraction  steps  have  a  o(l)  probability  of 
snccess  only  after  the  component  size  has  been  reduced  to 
0(7i)  where  7i  is  the  genus  of  component  i. 

Theorem  3  After  O(iblgTi)  parallel  paired  contraction 
steps, any  connected,  degree-40,  graph  of  genus  7  has  con¬ 
tracted  to  a  graph  with  0{y)  vertices  with  probability  1  — 
0(l/n*)  for  any  constant  k. 

Proof.  During  each  parallel  contraction  step,  each  vertex 
chooses  an  edge  randomly  out  of  all  adjacent  edges  eligible 
for  contraction.  Since  in  phase  1  no  vertex  degree  exceeds 
40  throughout  the  contraction,  on  each  iteration  every  good 
vertex  «  has  at  least  a  1/40  probability  of  merging,  since 
1/40  is  a  lower  bound  on  the  probability  that  the  edge  («,  v) 
that  vertex  «  chooses  is  also  chosen  by  the  vertex  v  at  the 
other  end.  Lemma  2  shows  that  a  bounded-degree  genus- 
7  graph  of  size  4O7  will  have  at  least  1/12  of  its  vertices 
be  good  during  each  iteration  of  phase  1.  Therefore,  we 
expect  that  at  least  n/480  of  the  vertices  of  an  n-vertex 
bounded-degree  graph  will  contract  on  each  iteration  of 
phase  1  provided  that  n  is  sufficiently  large  compared  to  the 
genus  of  the  graph.  A  Chernoff-bound  argument  completes 
the  proof.  B 

We  have  shown  that  with  high  probability,  we  require 
only  O(lgn)  parallel  contraction  steps  to  reduce  connected 
component  t  to  size  0(ti)  and  hence  reduce  the  maxi¬ 
mum  connected  component  to  size  0(7).  Therefore  phase 
1  runs  in  <?(lgn)  time  and  with  high  probability  reduces 
the  largest  component  of  the  graph  to  size  0(7). 

Let  us  now  consider  the  graph  at  the  start  of  phase  2. 
Assuming  that  we  are  not  already  done,  we  have  with  high 
probability  at  most  0(7)  nodes  in  each  component.  Con¬ 
nected  subgraphs  contract  independently  so  the  asymptotic 
contraction  time  during  phase  2  is  dominated  by  the  time 
to  contract  the  largest  component.  If  the  maximum  com¬ 
ponent  has  size  s,  then  by  the  analysis  in  section  2,  phase  2 
'.erminates  in  O(lg’s)  time.  Since  we  have  s  =  0(7)  with 
:.igh  probability,  then  phase  2  terminates  in  0(lg^  7)  time 
with  high  probability. 

Since  the  time  to  perform  the  contraction  is  simply  the 
sum  of  the  times  to  perform  phase  1  and  phase  2,  we 
have  that  an  n-vertex  genus  7  graph  is  contracted  by  al¬ 
gorithm  Bounded-Contract  in  0(lg  n  -f-  Ig^  7)  time  with 
high  probability.  The  probability  of  algorithm  BoUNDED- 
CoNTRACT  failing  is  the  probability  of  failing  in  phase  1, 
since  phase  2  is  deterministic. 

The  deterministic  version  guarantees  that  phase  1  re¬ 
duces  the  n-vertex  graph  to  an  0(7)-vertex  graph,  where 
7  is  the  largest  genus  of  any  connnected  component. 
Lemma  2  guarantees  that  in  each  parallel  contraction  step, 
each  ni- vertex  genus  7,  component  for  which  n,  =  fl(7i) 
has  at  least  n,/12  good  vertices.  Of  these  good  vertices, 
at  least  2/41  are  matched  by  the  symmetry- breaking  tech¬ 
niques  presented  in  section  3.  The  worst  case  is  achieved 
by  many  vertices  with  40  degree-one  neighbors.  There¬ 
fore,  each  such  component  is  reduced  by  at  least  a  factor 


of  491/492  after  each  parallel  contraction  step.  By  running 
phase  1  of  algorithm  Bounded- Contract  for  lg4«3/49]  n 
iterations,  we  are  guaranteed  that  each  component  is  of  size 
0(^i)  when  we  proceed  to  phase  2.  In  phase  1,  the  time 
required  to  perform  each  parallel  contraction  step  is  dom¬ 
inated  by  the  @(lg*  n)  tim^  required  to  color  rooted  trees 
[11].  Therefore  phase  1  always  terminates  in  6(lgnlg*n) 
time.  Because  the  graph  passed  to  phase  2  is  always  of 
size  0(7),  phase  2  always  terminates  in  time  0(lg^  t).  and 
therefore  the  deterministic  contraction  algorithm  runs  in 
time  0(lg  tt  Ig*  n  -Hg*  7). 

There  is  a  constaat-iactor  tradeoff  between  time  spent 
in  phase  1  and  time  spent  in  phase  2.  The  constant  fac¬ 
tor  1/12  of  good  nodes  guaranteed  by  lemma  2  can  be  re¬ 
placed  by  1/c  for  any  constant  c  >  6.  The  more  general 
version  of  lemma  2  states  that  any  bounded-degree  genns-y 
graph  with  at  least  20c7/(c— 6)  nodes  has  at  least  1/c  good 
nodes.  If  we  choose  to  base  our  algorithm  upon  a  fraction 
of  1/c  >  1/12  good  nodes,  the  constant  j  in  line  1  of  algo¬ 
rithm  Bounded- Contract  decreases.  The  expected  size 
of  the  graph  passed  on  to  phase  2  increases,  however,  and 
therefore  we  can  expect  phase  2  to  require  m'>re  contraction 
steps. 

We  should  comment  that  although  the  constants  are 
large  in  the  asymptotic  bounds  in  Theorem  3,  the  anal¬ 
ysis  is  highly  pessimistic.  Typically,  a  vertex  has  a  much 
greater  chance  than  1  in  12  of  being  good,  and  if  it  is  good, 
it  typically  has  more  than  a  1  in  40  chance  of  merging  with  a 
neighbor,  because  the  neighbor  is  unlikely  to  be  incident  on 
40  eligible  edges.  Consequently,  in  practice  we  could  reduce 
the  constant  j  in  line  1  of  algorithm  Bounded- Contract 
without  significantly  harming  the  behavior  of  the  contrac¬ 
tion  algorithm. 


5  Missing-Edge  Lemma 

In  this  section  we  present  the  proof  of  the  missing-edge 
lemma  used  in  the  analysis  of  algorithm  Bounded  - 
Contract.  We  begin  by  defining  cycle  splitting  of  a  con¬ 
nected  genus-T  graph,  a  technique  that  will  be  used  in  the 
inductive  proof  of  the  missing-edge  lemma.  We  then  prove 
that  performing  cycle  splitting  on  a  genus-7  graph  results 
in  a  new  connected  graph  of  genus  strictly  less  than  7  or  re¬ 
sults  in  two  disjoint  graphs  whose  geni  sum  to  the  original 
genus  7.  Finally,  the  proof  of  the  missing-edge  lemma  is  a 
double  induction  on  genus  and  number  of  bad  degree-three 
vertices.  Throughout  this  section,  we  assume  all  graphs 
are  connected.  A  disconnected  graph  G  whose  genns  7a  is 
the  sum  of  the  geni  of  its  connect^  components  can  only 
be  sparser  than  a  connected  graph  of  genns  7a  and  hence 
the  disconnnected  graph  cannot  have  a  weaker  missing-edge 
lemma. 

Suppose  we  have  a  connected  graph  of  genns  7  embedded 
on  a  surface  of  genns  7  (e.g.  a  sphere  with  7  handles).  Con¬ 
sider  any  simple  cycle  (og ,  vi , . . .  V(_] )  of  the  graph.  As  we 
travel  along  the  cycle  from  vertex  tag  back  to  vertex  vg,  we 
have  a  well-defined  notion  of  right  and  left,  since  the  surface 
upon  which  the  graph  is  embedded  is  orientable  [7].  Thus 
the  vertices  adjacent  to  each  vertex  Vi  on  the  cycle  can  be 
partitioned  into  two  sets:  those  that  connect  to  the  vertex 
from  the  right  (Vi,n)  and  those  that  connect  to  the  vertex 
from  the  left  (V'i.r.).  We  perform  a  cycle  splitting  opera¬ 
tion  as  follows.  We  remove  the  cycle  from  the  graph,  split 
each  vertex  on  the  cycle  into  two  vertices  v,,k  and  v,  i 
and  form  them  into  two  cycles  (vo,A,vi,n . and 


Figure  3;  Given  an  embedding  of  a  graph,  a  cycle  can 
be  split  into  two  pieces:  the  right  piece  which  is  con¬ 
nected  to  nodes  on  the  right  as  we  tranvene  the  cycle, 
and  the  left  piece  which  is  connected  to  nodes  on  the 
left. 


(vo.i.tij.t, . .  ■  .vi-i.t).  Finally,  vertex  Vi,R  is  connected  to 
ea^  vertex  in  Vt.jt,  and  vertex  «i,x,  is  connected  to  each 
vertex  in  Vi,z,.  Intuitively,  if  we  view  the  embedded  cycle 
as  having  finite  thickness,  we  simply  cut  it  in  half.  In  fig¬ 
ure  3,  cycle  «o, . .  ■  «3  is  split.  The  left  version  is  connected 
to  edges  entering  the  cyde  from  the  left  and  vice  versa. 

The  following  lemma  wiU  be  used  to  allow  apidiution  of 
the  induction  hypothesis  in  the  proof  of  the  missiag-edge 
lemma. 

Lemma  4  //  we  split  o  cycle  C  s  of  o 

graph  G  of  genu*  7  to  obtain  a  new  graph  G'  of  gemu  7', 
we  have  either 

1.  graph  G'  contains  two  disjoint  eompoapnts:  graph  Gj, 
of  penus  7(,  and  graph  Gr  of  penus  in  such  that  71  + 
7s  *  7.  or 

i.  graph  G'  is  connected  and  7' s  7  —  1 . 

Proof.  Suppose  we  have  an  embedding  of  graph  G  os  an 
orientable  manifold  of  genus  7  such  as  a  sphere  with  7 
handles.  By  definition  of  the  genus  of  a  graph  and  of  a 
surface,  all  the  faces  of  graph  G  are  simply  connected  [T]. 
That  is,  they  can  be  continuously  contract^  to  a  point  on 
the  surface.  Let  v,  e,  and  /  be  the  number  of  vertices, 
edges,  and  faces  respectively  of  graph  G.  Then  the  Euler 
Characteristic  of  the  surface  is  defined  asx=v  —  e  +  /.  If 


Figure  4:  Given  an  embedding  of  a  graph  G  cm  a  sur¬ 
face  of  corresponding  genus,  we  can  remove  a  simple 
cycle  of  graph  G  from  the  surface  and  patch  the  re¬ 
sulting  hcdes  with  disks.  The  result  is  an  embedding 
of  the  graph  G'  obtained  by  splitting  the  given  cycle 
of  G.  The  new  surface  may  or  may  not  be  connected. 


we  embed  any  graph  on  the  surface  such  that  all  faces  are 
simply  connect^,  the  quantity  «  —  e  +  /  for  that  graph  will 
always  be  equal  to  the  Euler  characteristic  of  the  surface. 
Furthermore,  by  definition,  we  have  that  x  =  2  — 27  where 
7  is  the  genus  of  the  surface  (and  of  any  graph  which  can 
be  embedded  on  that  surface  such  that  each  face  is  simply 
connected). 

Consider  now  the  surface  with  the  embedding  of  graph 
G.  We  cut  the  surface  along  cycle  C  and  patch  the  two 
resulting  boundaries  with  disks  as  illustrate  in  figure  4. 
The  resulting  surface  is  an  orientaUe  manifold  that  may  or 
may  not  be  connected. 

Let  ns  caknlate  the  Euler  characteristic  of  the  new  snr- 
fooe.  The  new  graph  G'  vs  embedded  on  the  new  surface 
such  that  all  faces  are  simply  connected  since  the  disks  used 
to  patch  the  cut  are  simply  connected  and  no  other  faces  of 
the  original  graph  G  are  ^tered  by  the  procedure.  There¬ 
fore,  if  o',  e',  and  /'  are  the  number  of  vertices,  edges,  and 
faces  in  the  cycle-split  graph  G',  then  the  Euler  characteris¬ 
tic  of  the  new  snrf^e  is  equal  to  x'  v'-e'-f /'  =  2-27' 
where  7'  is  the  genns  of  graph  (j  and  the  surface  upon 
which  it  is  embedded.  We  have  that  v' s  o-fl  and  e' s  e-f  ( 
and  /'  s  /  -t-  2.  Therefore,  we  have  that  x'  —  X  +  3. 

First  we  will  consider  the  case  where  the  new  surface  is 
disconnected.  In  this  case  the  Euler  characteristic  is  sim|dy 
the  sum  of  the  Enkt  characteristics  of  the  two  pieces  (the 
formula  v  —  e  -I-  /  is  additive).  Let  one  piece  have  character¬ 
istic  Xr  =  2  —  27 1  and  the  other  piece  have  characteristic 
Xn  s  2  -  27n.  Tben  we  have  that 

X'  »  Xn  +  XL 

=s  2  -  27i  +  2  -  27h 
s=  4  -  2(71,  7«). 

We  also  have  that  x'  *  X  +  2  s  4  —  27.  Therefore,  we  have 
that 

4-27  =  4  -  2(71.  +  7s) 

7  =  71.  +  7s 

which  proves  case  (f). 

Next  we  consider  the  case  where  the  new  surface  is  con¬ 
nected.  Then  we  have  that  x'  =  2  —  27'  and  x'  =  X  +  2  = 
4  —  27-  Therefore,  equating  the  right-hand  sides,  we  have 


thAt  2  -  27' s  4  —  27  or  7'  =  7  —  1  wkich  prove*  caae  (2). 

■ 

Now  we  proceed  to  the  proof  of  the  nuanng-edge  lemma. 
First,  let  ns  remind  the  reader  of  tome  terminology.  Bad 
vertices  are  defined  as  they  were  in  section  4:  no  adjacent 
edge  can  be  contracted  without  potentially  creating  a  ver¬ 
tex  with  degree  exceeding  the  bound.  Bad  degree-three 
vertices  must  be  adjacent  to  three  degree-dmu  indepen¬ 
dent  vertices  (no  edges  between  any  pair).  Thu*  degree- 
three  vertke*  force  a  stronger  sparsity  than  Euler’s  for¬ 
mula  alone  since  these  three  muting  edges  can  be  added  to 
the  graph  without  increasing  its  genus.  The  missiag-edge 
lemma  quantifies  this  increased  sparsity. 

For  convenience  we  restate  the  lemma  as  it  appeared 
in  section  4,  assuming  this  time  that  the  graph  is  con¬ 
nected.  The  proof  goes  through  for  disconnected  graph* 
where  genus  is  defied  in  the  usual  way. 

Leminn  1  (Missing-Edge  Lemma)  A  graph  G  =  (V,  £} 
of  genu*  7  with  63  >  0  bad  degree-3  vertiee*  ha*  at  mo*t 
3|V|  —  63  -f  IO7  edge*. 

Proof.  The  proof  centers  on  showing  that  the  number  of 
missing  edges  is  at  least  63  —  47  -f  1  for  63  >  0,  which, 
together  with  the  constraint  that  |i^|  <  3  |V|  -i-  67  from 
Elder’s  formula,  suffices  to  prove  the  lemma.  The  goal  of 
the  proof  is  to  show  that  sharing  of  missing  edge*  is  limited. 
For  the  remainder  of  this  proof,  when  we  say  “bad  vertex,” 
we  assume  the  vertex  has  degree  3. 

We  begin  our  induction  by  showing  that  the  lemma  holds 
for  the  caae  7  s  0  for  any  number  63  of  bad  vertices.  That 
is,  we  show  that  a  planar  graph  G  =  (V,  E)  with  63  >  0 
bad  degree-3  vertices  has  at  most  3 1  V|  —  ba  edges. 

The  proof  of  the  planar  case  is  by  induction  on  ba.  If 
^3  —  0,  Euler’s  formula  alone  is  sufficient.  The  lemma 
holds  trivially  for  the  case*  is  =  1  and  is  —  2  since  any 
graph  with  at  least  one  bad  degree-3  vertex  has  at  least 
three  distinct  missing  edges.  Now  consider  the  case  is  =  3. 
if  there  are  only  three  missing  edge*  in  the  graph,  then  each 
of  the  three  baid  vertices  is  associated  with  the  same  three 
missing  edge*.  Hence  all  three  bad  vertices  are  adjacent  to 
the  same  three  degree-dm*!  vertices,  and  hence  the  graph 
contains  an  instance  of  Ka.a,  which  violate*  the  assumption 
that  the  graph  is  planar.  Consequently,  the  graph  has  at 
least  4  =  is  -f- 1  missing  edge*. 

We  now  consider  the  general  planar  case  is  >  4.  Assume 
inductively  that  any  graph  with  i  <  is  bad  vertices  has  at 
least  k  +  i  missing  edges,  bat  that  there  exists  a  planar 
graph  G  w  (V,  E)  with  is  bad  vertices  and  m  <  is  missing 
edges.  Since  each  bad  vertex  is  associated  with  exactly 
three  missing  edges,  it  follows  that  there  exist*  a  missing 
edge  («i,«3)  associated  with  at  least  three  bad  vertices 
V],  «3,  and  «3.  There  must  be  at  least  one  additional  bad 
vertex  ««,  since  we  have  is  >  4.  For  some  embedding  of  G 
on  the  sphere,  the  Jordan  Cnrve  Theorem  ensures  that  two 
of  the  vertices,  say  «i  and  vj,  together  with  m  and  ns  form 
a  cycle  C  that  separates  from  V4,  as  shown  in  Figure  S. 

We  are  now  set  up  to  apply  the  induction  hypothesis. 
Let  Gim  be  the  subgraph  of  G  induced  by  cycle  C  and  the 
vertices  on  one  side  of  C,  and  let  Gn\  be  the  subgraph 
induced  by  cycle  C  and  the  vertices  on  the  other  side  of 
C.  Add  to  m  and  ns  in  each  of  Gi.  and  G^.t  enough 
degree-I  neighbor*  to  maintain  their  degrees  a*  dm**-  Let 
btm  be  the  number  of  bad  (degree-3)  vertices  in  Gii ,  and  let 
i»at  be  the  number  of  bad  vertices  in  Gout-  Since  no  new 
bad  degree-3  vertices  are  introduced  into  Gu  and  G«at,  we 
have  6is  <  ^3  and  6o«i  <63-  By  the  induction  hypothesis, 
therefore,  graph  Gm  has  at  least  is.  -f- 1  missing  edges,  and 


Figure  5:  A  cycle  C  =  (ui,  vi, us,  vs)  that  separates  G 
into  two  subgraphs  Gin  uid  Gout>  ^th  at  most 
is  —  1  bad  degree-3  vertices. 


graph  Gmi  has  at  least  b,mt  +  1  musing  edge*. 

We  next  account  for  interaction*  between  graphs  Gi.  and 
Goat  to  obtain  the  lower  bound  of  is  +  1  on  the  number  of 
missiag  edge*  in  the  original  graph  G.  First,  observe  that 
hia  +  bent  =  ba,  since  when  oi  is  a  bad  vertex  in  one  of  the 
graphs,  it  is  a  degree-2  vertex  in  the  other,  and  similarly 
for  02-  Moreover,  each  of  oi  and  vj  are  Lad  in  at  least  one 
of  Gin  and  Goat  becanse  we  dummied  up  each  of  sj  and  *} 
to  have  degree  dmai-  The  number  of  irissing  edges  in  G  is 
at  least  (bn  +  1)  -f  (bent  +  1)  minus  the  number  of  missing 
edges  shared  by  Gi,  and  Goat-  The  number  of  such  shared 
missiag  edges  is  in  fact  1,  namely,  the  edge  (si, *3),  since 
«i  and  *3  are  the  only  degree-daiu  vertices  in  both  Gia  and 
Gent-  Consequently,  the  number  of  missing  edges  in  G  is 
at  least  (lia  4-1)4-  (ioai  4- 1)  —  1  =  is  -f  1,  which  complete* 
proof  of  the  ptaaar  version  of  the  lemma. 

We  have  shown  that  the  missing  edge  lemma  lemma 
holds  for  the  case  7  »  0  for  any  number  ba  of  bad  ver¬ 
tices.  To  complete  the  base  cases,  we  see  that  for  genera] 
graphs  the  lemma  holds  for  is  —  0  trivially  using  only  Eu¬ 
ler’s  formula.  It  also  holds  trivially  for  1  <  is  <  47  4-  3  for 
all  7  since  then  the  lemma  only  require*  that  the  number 
of  missing  edges  is  at  least  3  which  is  true  for  any  graph 
with  at  least  one  bad  vertex. 

We  now  consider  the  general  case  is  >47-1-3  and  7  >  1  ■ 
We  assume  inductively  that  any  genns-7  graph  with  k  <  ba 
bad  vertices  has  at  least  i  —  47  -t- 1  missing  edges  and  any 
graph  with  genus  >  <  7  and  any  number  is  bad  degree- 
three  nodes  has  at  least  is  —  4y  -f  1  missing  edges.  Assume, 
however  that  there  exists  a  graph  G  s  (V,  E)  of  genns  7 
with  is  >47-1-4  bad  vertices  and  m  <  63  —  47  missing 
edge*.  Since  each  bad  vertex  is  associated  with  exactly 
three  edges,  it  fellow*  that  there  exists  a  missing 

edge  (m,**)  associated  with  at  least  three  bad  vertices 
«!,  V3,  and  «3.  There  mnst  be  at  least  one  additional  bad 
vertex  ««  since  we  have  ba  >7.  The  situation  is  illustrated 
in  figure  5. 

To  allow  application  of  the  induction  hypothesis,  we  split 
graph  G  along  one  of  the  three  simple  cydes  shown  in  fig¬ 
ure  5  to  obtain  a  new  graph  G'  of  genns  7'.  As  we  did  in 
proving  the  planar  version  of  this  lemma,  we  add  degree- 
one  neighbors  to  the  high-degree  nodes  «i  and  *3  so  that 
each  of  the  four  node*  resulting  from  the  split  of  si  and  *3 
has  maximum  degree.  Each  bad  degree-three  node  on  the 
cyde  is  split  into  two  nodes,  but  exactly  one  of  these  two 
nodes  is  a  bad  degree-three  node  in  the  new  graph  G'.  The 
other  node  split  from  that  vertex  is  of  degree  2.  Therefore 
graph  G',  whether  connected  or  not,  has  exactly  63  bad 


node*. 

Suppose  we  cut  split  grsph  G  along  oon.  of  the  cycles 
such  that  the  resulting  graph  G'  is  stiU  connect^.  Then 
by  lemma  4(f)  we  have  that  7'  =  7  —  1.  As  argned  above, 
graph  G'  has  the  same  number  63  of  bad  degree-three  ver¬ 
tices  as  the  original  graph  did.  Because  graph  G'  has  lower 
genus,  we  can  apply  the  induction  hypothesis.  Therefore, 
graph  G'  has  at  least  63  —  4(7  —  1)  -f  1  =  6s  —  47  -H  5 
missing  edges.  We  overconnted  exactly  one  missing  edge, 
namely  («i,ii3),  since  these  are  the  only  two  degree-dm^i 
vertices  duplicated.  Therefore  we  have  at  least  63  -47-1-4  > 
63—47-1-1  missing  edges. 

Now  suppose  that  splitting  along  any  of  the  three  simple 
cycles  separates  G'  into  two  graphs:  graph  Gl  with  genus 
7l  and  bi  bad  degree-three  vertices,  and  graph  Gn  with 
genus  7/t  and  6n  bad  degree-three  vertices.  Suppose  that 
splitting  along  one  of  the  cycles  yields  graphs  such  that 
71,  <  7  and  7n  <  7.  Applying  the  induction  hypothesis  to 
each  side  we  have  that  the  number  of  missing  edges  in  graph 
Gr  is  at  least  hn— 47n-l-l  and  the  number  of  missing  edges 
in  graph  Gi  is  at  least  bi  —  ^"II  1.  Adding  the  missing 
edges  together  and  subtracting  one  for  the  overconnt  of 
edge  (uj ,  ttj),  we  have  that  the  original  graph  G  has  at  least 
(6L-l-6ii)-4(7i.-l-7ii)+2-l  missing  edges.  From  lemma  4(/) 
we  have  that  71  -1-  7/t  =  7  and  we  argued  earlier  that  6x,  -1- 
hn  =  63 .  Substituting  these  back  into  the  expression  for  the 
number  of  missing  edges  in  graph  G,  we  have  that  graph 
G  has  at  least  63  —  47  -1- 1  missing  edges. 

The  final  case  we  must  consider  is  that  each  of  the  three 
cycles  cuts  off  a  planar  patch.  In  this  case  we  have  a  sit¬ 
uation  analogous  to  the  planar  case  argued  earlier.  As  il¬ 
lustrated  in  figure  5,  one  of  the  cycles  cuts  the  graph  into 
two  graphs  Gi  and  Gr  such  that  6x,  <  63  and  bR  <  63. 
Without  loss  of  generality,  let  us  assume  that  graph  Gr  is 
planar.  Then  by  lemma  4(1),  graph  Gl  must  have  genus 
7-  Applying  the  induction  hypothesis  to  the  two  graphs 
(since  each  has  less  than  63  bad  vertices),  we  have  that 
Gr  has  at  least  6r  -f  1  missing  edges  and  Gl  has  at  least 
6l  -47‘f  1  missing  edges.  Again  we  have  that  bR-^bL  ~  bs. 
Adding  the  missing  edges  and  subtracting  one  for  the  over¬ 
count  of  edge  (ui,tt3),  we  have  that  graph  G  has  at  least 
6l  -f  6n  —  47  -f  2  -  1  =  63  -  47  -f  1  missing  edges  which 
completes  the  proof.  B 


6  Applications 

The  most  direct  application  of  the  parallel  contraction  algo¬ 
rithm  is  to  the  vision  problem  of  determining  the  connected 
regions  of  an  image  represented  as  a  two-oimensional  ar¬ 
ray  of  pixels.  We  can  view  the  image  as  a  planar  graph 
G  s  (V,  £),  where  the  vertex  set  V  is  the  set  of  pixels,  and 
(a,v)  €  if  *  and  «  are  adjacent  pixels  with  the  same 
color.  The  region  labeling  problem  is  to  assign  each  pixel 
an  integer  such  that  two  pixels  have  the  same  label  if  and 
only  if  they  are  path  connected  in  G. 

Region  labeling  can  be  solved  by  a  connected- 
components  algorithm.  Shiloach  and  Vishkin  [27]  give  an 
0(lgn)-time,  n-processor  parallel  algorithm  for  connected 
components  using  the  concurrent- read,  concurrent-write 
(CRCW)  PRAM  model.  Hagerup  [13]  gives  an  O(lgn)- 
time,  n/lg  n-processor  CRCW  algorithm  for  graphs  of 
bounded  genus*.  The  best  algorithm  to  date  in  the 

’  HagCTup  has  recently  independently  devel¬ 
oped  an  0(lgnlg*  n)-time  n/Ig nig*  n-processor  deterministic 
EREW  algorithm  for  graphs  of  bounded  genus  [14] 


more  easily  implemented  exclusive-read,  exclusive-write 
(EREW)  model  is  due  to  Hirschberg,  Chandra,  and  Sar- 
wate  [17],  who  give  an  G(lg^  n)  time  algorithm  for  con¬ 
nected  components  using  the  adjacency-matrix  representa¬ 
tion  of  a  graph.  Leiserson  and  Maggs  [20]  give  an  0(lg^  n) 
step,  n-processor,  randomised  connect^-components  algo¬ 
rithm  for  a  DRAM  (distributed  random-access  machine), 
an  EREW -like  model  that  includes  the  cost  of  communica¬ 
tion.  Blefloch  [2]  gives  an  G(lg  n)  randomised  algorithm  for 
a  model  that  includes  parallel  prefix  as  a  primitive  opera¬ 
tion.  Lim  [21]  gives  a  region-labidiag  algorithm  that  runs  in 
0(lg^  n)  time  on  an  EREW  PRAM  which  uses  the  planar 
and  geometric  properties  of  a  mesh.  Gasit  and  Reif  have 
recently  developed  a  randomised  0((n-f  m)/  Ig  n)-procesaot 
EREW  algorithm  that  runs  in  time  0(lg  n  -f  Ig^  g)  where 
g  is  the  genus  of  the  graph  (what  we  call  70,  which  can  be 
much  larger  than  7  in  disconnected  graphs).  They  require 
an  embedding  of  the  graph.  Other  algorithms  for  connected 
components  are  given  by  [6,  9,  19,  24,  30,  32]. 

Our  algorithm  for  labeling  planar  graphs  is  asymp¬ 
totically  the  fastest  algorithm  to  date  in  the  EREW 
PRAM  model.  The  algorithm  uses  Algorithm  Bounded- 
Contract  to  simultaneously  contract  each  component  of 
the  graph  to  a  single  node,  and  then  it  simply  reverses 
the  contraction  process  and  assigns  the  same  label  to  all 
nodes  in  each  connected  component.  The  algorithm  uses 
0(n/ lgn)processors,  and  with  high  probability  (i.e.,  prob¬ 
ability  at  least  1  —  0(l/n*)  for  any  constant  k),  the  al¬ 
gorithm  runs  in  G(lgn)  time  on  an  EREW  PRAM  or  a 
DRAM.  Thus,  onr  algorithm  achieves  the  best  possible 
asymptotic  bound  in  the  most  restrictive  parallel  mod¬ 
els.  The  deterministic  version  of  the  algorithm  runs  in 
0(lg  n  Ig*  n)  time  using  an  optimal  number  of  processors 
(n/lg nig*  n).  Like  Lim’s  0(lg’  n)-time  algorithm,  our  al¬ 
gorithm  for  region  labeling  takes  advantage  of  the  planar 
nature  of  the  image  graph,  but  unlike  his  algorithm,  it  does 
not  depend  on  the  geometric  nature  of  the  image.’  Thus, 
our  algorithm  works  on  irregular  planar  structures,  includ¬ 
ing  meshes  with  local  refinements  and  meshes  with  (static) 
faulty  elements. 

The  contraction  algorithm  can  also  be  used  for  a  vari¬ 
ety  of  other  applications  including  algorithms  for  bicon- 
nected  components  of  planar  graphs  and  for  spanning  trees 
of  planar  paphs.  The  running  times  for  these  algorithnu  is 
asymptotically  the  same  as  for  the  contraction  algorithm, 
and  for  these  problems  the  running  times  are  asymptoti¬ 
cally  the  best  to  date  in  the  EREW  model.  All  the  al¬ 
gorithms,  including  the  connected  components  algorithm, 
can  be  generalised  to  paphs  of  higher  genus  and  higher 
degree  using  algorithm  Bounded-Contract.  The  asymp¬ 
totic  running  times  for  these  problems  match  the  running 
times  of  the  contraction  algorithm  (i.  e.  0(lg  n-f-lg’  7)  with 
high  probability  where  7  is  the  maximum  genus  of  any  con¬ 
nected  component).  For  classes  of  graphs  with  genus  o(n) 
this  is  the  best  bound  to  date.  For  graphs_^of  genus  n(n) 
this  matches  the  best  previous  bound  of  0(lg’  n). 


’Subsequent  to  the  development  of  the  planar  version  of 
our  algorithm,  Lim,  Apawal,  and  Nekludova  [32]  obtained  an 
0(lgn)-time  deterministic  algorithm  for  region  labeling.  Like 
Lim’s  earlier  work,  their  algorithm  exploits  the  geometric  prop¬ 
erties  of  specific  plattar  grids. 


7  Conclusion 

Ib  this  Mction  we  comment  npon  the  pn^cnl'behnTior 
of  algorithm  BoUNDElvComnACTand  preaent  lorae  open 
probiema. 

In  aection  4  we  commented  that  the  behavior  of  algo¬ 
rithm  Bounded- CoNT*'^CTin  practice  ia  likely  to  be  mnch 
better  than  that  gnarenteed  by  our  worst-caae  analyaia. 
Thia  ia  particolaxlv  true  for  the  bounded-genua  caae  where 
we  explicitly  r*  zk  for  termination.  For  example,  we  arc 
not  likely  to  require  the  worat-caae  number  of  iterationa 
to  reduce  each  component  of  a  planar  graph  to  a  aingle 
node.  The  general  algorithm,  however,  alwaya  performa 
the  worat-caae  number  of  iterationa  for  phaae  1  becauae  we 
do  not  have  a  meana  of  detecting  when  phaae  1  ia  com¬ 
pleted  nnleaa  we  in  fact  complete  the  entire  contraction. 
Phaae  2  only  performa  aa  many  contraction  atepa  aa  necea- 
aary.  Becauae  of  thia  control  atmctnre,  it  is  probably  good 
in  practice  to  chooae  the  constant  j  in  line  1  of  algorithm 
Contract  baaed  upon  lemma  2  with  a  fraction  of  1/c  good 
nodes  for  constant  c  =  6  -f  <  and  <  >  0. 

To  detect  the  termination  of  phaae  1  we  require  an 
0(lgn)  time  n  processor  algorithm  which  can  determine 
the  genua  of  an  n  node  graph  to  within  a  constant  with¬ 
out  knowing  the  value  of  n.  The  problem  appears  difficult 
since  determining  the  exact  genua  7  of  an  n-node  graph  is 
NP-complete  [29]  and  the  beat  sequential  algorithm  for  this 
problem  runs  in  time  [8].  There  aeema  to  be  no 

work  on  parallel  algorithms  which  approximate  the  genus 
of  a  graph. 
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